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■ Abstract 

We start with a set of n players. With some probabihty P{n, k), we kiU n — k players; 
the other ones stay alive, and we repeat with them. What is the distribution of the number 
Xn of phases (or rounds) before getting only one player? We present a probabilistic 
analysis of this algorithm under some conditions on the probability distributions P{n, fc), 
including stochastic monotonicity and the assumption that roughly a fixed proportion a 
of the players survive in each round. 
t/3 , We prove a kind of convergence in distribution for X„ — \ogi/^ n; as in many other 
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similar problems there are oscillations and no true limit distribution, but suitable sub- 
sequences converge, and there is an absolutely continuous random variable Z such that 
d{Xn, \Z+\ogj^/^ n] ) — > 0, where d is either the total variation distance or the Wasserstein 
^ , distance. 

Applications of the general result include the leader election algorithm where players 
, are eliminated by independent coin tosses and a variation of the leader election algorithm 

proposed by W.R. Franklin [7]. We study the latter algorithm further, including numerical 
, results. 

o 

oo . 

O ■ 1 A general convergence theorem 



We consider a general leader election algorithm of the following type: We are given some 
random procedure that, given any set of n > 2 individuals, eliminates some (but not ah) 
individuals. If there is more that one survivor, we repeat the procedure with the set of 
survivors until only one (the winner) remains. We are interested in the (random) number Xn 
of rounds required if we start with n individuals. (We set Xi = 0, and have Xn > 1 for n > 2.) 
We let N^. be the number of individuals remaining after round k; thus Xn '■= min{/c : Nf^ = 1}, 
where we start with A'^o = n. For convenience we may suppose that we continue with infinitely 
many rounds where nothing happens; thus A^ is defined for all > and N/. = 1 for all 
k > Xn. 

We assume that the number Yn of survivors of a set of n individuals has a distribution 
depending only on n. We have 1 < 1^ < n; we allow the possibility that Yn = n, but we 
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assume P(l^ = n) < 1 for every n > 2, so that we will not get stuck before selecting a 
winner. We further assume that, given the number of remaining individuals at the start 
of a new round, the number of survivors is independent of the previous history. In other 
words, the sequence (Nf^)^ is a Markov chain on {1, 2, . . .}, and X„ is the number of steps 
to absorption in 1. The transition probabilities of this Markov chain are, with Yi = 1, 



P{i,j) := ¥{Yi = j) = P(j survives of a set of i) 



(1.1) 



Note that P{i,j) = if j > i and P{i,i) < 1, for i > 1. Conversely, any Markov chain on 
{1, 2, . . .} with such P{i,j) can be regarded as a leader election algorithm in the generality 
just described. 

We will in this paper treat leader election algorithms where, asymptotically, a fixed pro- 
portion is eliminated in each round. (Thus, we expect X„ to be of the order logn.) More 
precisely, we assume the following for Yn, where we also repeat the key assumptions above. 
(Here and below, logn should be interpreted as some fixed positive number when n = 1.) 

Condition 1.1. For every n > 1, Yn is a random variable such that 1 < 1^ < n, and 
P(y„ = n) <l for n>2. Further: 

(i) Yn is stochastically increasing in n, i.e., P(yn ^ k) > P(l^-|-i < k) for all n > 1 and 
k >1. Equivalently, we may couple Yn and Yn+i such that Yn < Yn+i- 

(ii) For some constants a G (0, 1) and e > and a sequence 6n = 0((log n)^^^^) , 

EYn+i-EYn = a + Oi5n). (1.2) 



(iii) For some e and 5n as in (ii) 

¥{\Yn - an\ > (5„n) = 0{n-^-'). 



(1.3) 



Note that 

for some p > 4 suffices for 
and e small). 



an ■ 



0(nP/2) 



m 



for a suitable choice of e > and (5„ (e.g., 6„ 



(1.4) 

n~^, rj > 



Remark 1.2. If (jl.2p or (jl.3p holds for some sequence ((5„), it holds for every larger sequence 
{6n) too; similarly, if 5n = 0((logn 



^ ^) or (II. Sp holds for some e, it holds for every smaller 
hold with the same £ > and the same 



m 



e too. Hence we may assume that (ii) and 

(5„, and we may assume 5n > (logn)~^~^. In particular, this implies that 6k = 0{6n) when 
C^^n < k < Cn, for each constant C. 



The behaviour of the election algorithm is given by the recursion Xi = and 



X„ 



Xy„ + 1, n > 2, 



(1.5) 



where we assume that (^i)"=i and Yn are independent. We state a general convergence 
theorem for leader election algorithms of this type. 
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We recall the definitions of the total variation distance dxv ^-nd the Wasserstein distance 
dw (also known as the Dudley, Fortet-Mourier or Kantorovich distance, or minimal Li dis- 
tance); these are both metrics on spaces of probability distributions, but it is convenient to 
write also d^y^X, Y) := d^y{n, u) and d-w{X, Y) := d-witJ-, i^) for random variables X, Y with 
X ~ /i and y ~ I/. 

The total variation distance c^tv between (the distributions of) two arbitrary random 
variables X and Y is defined by 

dTy{X, Y) := sup \F{X e A) -F{Y £ A)\. (1.6) 
A 

For integer-valued random variables, as is the case in our theorem, this is easily seen to be 
equivalent to 

dTv{X, Y) = IY1 = k)-FiY = k)\. (1.7) 

k 

Further, for any distributions /x and v, 

d^yifi, v) := inf {F{X ^Y) : X ^ ^l,Y ^ v] ; (1.8) 

the infimum is taken over all random vectors (X, Y) on a joint probability space with the 
given marginal distributions /i and v. (In other words, over all couplings {X, Y) ol ^ and 
v.) For integer- valued random variables, convergence in d^^y is equivalent to convergence in 
distribution, or equivalently, weak convergence of the corresponding distributions. 

The Wasserstein distance dw is defined only for probability distributions with finite ex- 
pectation, and can be defined by, in analogy with (|1.8p . 

dy^{^, u) ■= inf {E|X -Y\: X ^ ^^,Y ^ v] . (1.9) 

There are several equivalent formulas. For example, for integer-valued random variables, 

dTv{X,Y) = ^|P(X < k)-F{Y < k)\. (1.10) 
k 

It is immediate from (jl.Sp and (jl.9p that for integer-valued random variables X and Y (but 
not in general), 

dTY{X,Y)<dw{X,Y). (1.11) 

It is well-known that dw is a complete metric on the space of probability measures on M 
with finite expectation, and that convergence in dy/ is equivalent to weak convergence plus 
convergence of the first absolute moment. 

All unspecified limits in this paper are as n ^ oo. 

Theorem 1.3. Consider the leader election algorithm described above, with Yn satisfying 
Condition Then, there exists a distribution function F with bounded density function 

f = F' such that 

sup \F{Xn <k)- F{k - n)| ^ (1.12) 

or, equivalently, if Z ^ F, 

dTv(^n, r^ + logi/„nl) ^0. (1.13) 
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More precisely, d-w{Xn, \Z + logj^/Q, n] ) ^ 0, which is equivalent to 

Y,\nXn<k)-Fik-logy^n)\^0. 



As a consequence, defining A.F{x) := F{x) — F{x — 1), 

sup \¥{Xn = k)- AF{k - logua n)\ 

Furthermore, 



EXr, 



:i.i4) 



(1.15) 



:i.i6) 



^ogi/an + (t>{n) + o{l), 

for a continuous function (j){t) on (0, oo) which is periodic in log]^/Q,t, i.e. (j){t) = (j){at), and 
locally Lipschitz. 

We thus do not have convergence in distribution as n — > oo, but the usual type of oscilla- 
tions with an asymptotic periodicity in log^/Q, n and convergence in distribution along subse- 
quences such that the fractional part {logj^/^^^ n} converges. (This phenomenon is well-known 
for many other problems with integer- valued random variables, see for example |141 llOj : it 
happens frequently when the variance stays bounded.) This is illustrated in Figure [H 




Figure 1: Illustration of Theorem 11.31 



Proof. We assume that (5„ are as in Remark II. 2 [ 
Let q := sup^>2P(l^n = n). Since each P(y„ 



n 



< 1, and P(y„ 



n) 



by (iii) 



q < 1. Hence Xn is stochastically dominated by a sum of n — 1 geometric Ge(l — q) random 
variables, and thus EX„ = 0{n). In particular, EX„ < oo for every n. 

Since the sequence (1^) is stochastically increasing, we may couple all such that 
Yi < Y2 < . . . ■ If we consider starting our algorithm with different initial values, and 
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use this coupling of (1^) in each round, we obtain a coupling of all X„, n > 1, such that 
^n+i > a.s. for every n > 1. We use these couplings of {Y^) and (X„) throughout the 
proof. 
Let 



dn 



EXn+1 - EXn = Xn+1 - Xr. 

max kdk- 

l<k<n 



We extend 6„ to real arguments by the same formula; thus, bt = for real t > 1. 
By dLl]), 

Xn = EXn = l+ExY„, n > 2. 



Thus, for n > 2, 



^n + l-l 



Y,djnYn<j <Yr,+i) 



By (ii) , KYn+i — Ey„ — > a, and thus there exists no such that if n > no then 



(1.17) 



Y^nYn < j < ^n+i) = Ey„+i - Ey„ < 1. 

Hence p.l7p implies, with d* = maxfc<n dfc, for n > no, 

dn{l - HYn <n< Yn+l)) < Y^^iYn <j< ^n+lX.i < <_i(l - F{Yn < n < 

i=i 

and thus < d^_^ so d*^ = dn\/ = c^n-i- Consequently, = (i*^^ < oo, for all n > no- 
In other words, d* := sup„ dn < oo. 

Let (3 := (1 + q)/2; thus q < /3 < 1. If n is large enough, so that {a + 5n+i){n + l) < (3n, 
then (fTTTl) yields, using ([L3]) and (p^l) . 



j<{a-S„)n ^ i j>/3n 



0{n 



+ - 



n(a - 5„) 
la + 0(5„)^ 



bpni^Yn+l-EYn) 



n a — 6n 



Thus 



ndn < (1 + 0(5n))&/3n + 0(n- 



(1.18) 
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Replace nhy k and take the supremum over all k such that /5n < k < n. Since bk is increasing, 
and by our simplifying assumptions in Remark 11.21 this yields 

6„ < (1 + 0{6n))bpn + 0{n-') = (1 + 0{6n))bpn- 

It follows by induction over m that if (l//?)™" < n < then 

m „ 

''"S<'.n(i+ig7) 

and thus 6„ = 0(1). In other words, we have shown 

dn = 0{l/n). (1.19) 

We now use the Wasserstein distance dw- Since Xn+i > Xn a.s., it is easily seen by (|1.9p 
that d^fil{Xn,Xn+l) = E{Xn+i - Xn) = dn- Thus, if m < n, by (|1.19p . 

n— 1 n— 1 / 

n — m 



dwiXn, X-ra) < ^ (-'^fc ; -''^fc+l ) = ^ dk = O 



k=m k=m 

and thus, for all n and m. 



dw(^n,X^) = 0(^V^)- (1-20) 

' n Am 



Note also that (iii) implies 

E\Yn - an\ < dnU + ©(n"^"^) = 0{n5n). (1.21) 

Define 

:=XLtj -logi/„t, t>l. (1.22) 
Then, for t > 2/a, using ([LS]), (fr20D . (fTni) . (fL3]) . and 1 < y^tj < 

(iw(-'^t,-'^Qt) = C^w(-'^[tJ - logl/a(i)i-'^M ~ logl/a(«*)) 

= o (t-iE|yLtj - M I) + o (tP (ntj < «V2)) 

= 0{6yt^ ) + 0{t- 1-^-') = O (log-i-= t) . 



Hence, for any t > 2/a, 



^ (iw(^„-n, ^a-.-ii) = O {log'' t) < oo. (1.23) 
j=0 
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Since dw is a complete metric, thus there exists for every t > a limiting distribution 
such that if Z{t) ~ then 

dv^{X^-,t,Z{t)) ^0 as j^oo. (1.24) 

In particular, 

X^-jt^Z{t) as j^oo. (1.25) 

(We find it more convenient to use the random variable Z{t) than its distribution /i(t).) 
Clearly, Z{at) = Z{t), so the distribution fi{t) is a periodic function of log^/^^j t. Hence, (|1.24p 
can also be written, adding the explicit estimate obtained from (ll.23p . 

dw {Xt, Z{t)) = O (log~^ t) ^ as t^oo. (1.26) 

Note further that, for 7 > 1, by ([02]) and (fOO]) . 

rfw(^t, < ^w(^LtJ , ^btj ) + I logl/a ^ - logl/a(7t)| = O ^^^^ ~ ^ + logi/, 7 

= 0(7-l + l/t). 

Replacing t by a~H and letting j — > cxd, it follows from (jl.24p that, for all t > and 7 > 1, 

d^{Z{t),Z{^t)) = 0{^-l). (1.27) 

Consequently, t — > ii{t) = C{Z{t)) is continuous and Lipschitz in the Wasserstein metric. 
Define, for every real x, 

F{x) = F{Z{t) < x) (1.28) 

for any t > such that x + logi/^ t is an integer; since Z{t) is periodic in log^/^ t, this does 
not depend on the choice of t. 

Since X^-jf- + logj^/^t = X^-^ — j G Z, the random variable Z{t) + logx/Q,i is integer- 
valued for every t. It is easily seen that for integer- valued random variables Zi and Z2, 
the total variation distance d'Y\i{Zi^ Z2) < dw(-2^i; ■^2)- Hence, for any x G M and u > 0, 
choosing t such that x + log^y^t is an integer and letting 7 = a~", which implies that 
X — u + logx/„(7t) = X + \ogii^{t) £ Z, we obtain from the definition (jl.28p and (|1.27p . 

F{x) - F{x -u)= F{Z{t) <x) - F(Z(7t) < x - u) 

= F{Z{t) + logi/„ t<x + logi/„ t) - P(Z(7t) + logi/„(7t) <x + logi/„ t) 

< dTY{Z{t) + logi/, t, Z{^t) + logi/„(7t)) 

< dw(^(t) + logi/„t, Z(7t) + logi/„(7t)) 

< (iw(^(i), ^(7*)) + I logi/, t - logi/„(7t)| 

= 0(7-l)+logi/„7 = 0(t^)- (1-29) 

Hence, F[x) is a continuous function of x. 

We have shown that t — > C{Z{t)) is continuous in the Wasserstein metric, and thus in 
the usual topology of weak convergence in the space 'P(R) of probability measures on M. 
Since further C{Z{t)) is periodic in t, the set {C{Z{t)) : t > 0} = {C{Z{t)) : 1 < t < a'^] 
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is compact in V(M), which by Prohorov's theorem means that the family {Z{t)} of random 
variables is tight, see e.g. Billingsley [1]. Hence, ¥{Z{t) < x) ^ as x ^ —oo and ¥(Z(t) < 
x) — > 1 as X ^ +00, uniformly in t, and it follows from (I1.28P that lim2._^_oo -^(a^) = and 
\imr,^ooF{x) = 1. 

Furthermore, ()1.26p and (11.280 show that, for any sequence kn of integers, as n ^ oo, 

< kn) = ¥{Xn < k^ - logi/„ n) = P(Z(n) <kn- logi/„ n) + o(l) 
= F(A;„ -logi/„ n) + o(l). 

Since further the sequence X„ is increasing, it now follows from Janson |10|, Lemma 4.6] that 
F is monotone, and thus a distribution function. By (|1.29p . the distribution is absolutely 
continuous and has a bounded density function F'{x). 

It is easy to see that ()1.30p . (|1.12p and (jl.lSp are equivalent, see flOl Lemma 4.1]. The 
corresponding result in the Wasserstein distance follows from (I1.26P because d-wiXn, \Z + 
logi/„n]) = dYj{Xn,Z{n)), e.g. by Remark [Ol below: (jl.l4p then follows by (jl.lOp . Finally, 
(jl.26p implies that 

EXt = EXl^j - logi/„ t = EZ{t) + O (log-^ t) , 

which proves (|1.16p with := EZ(t), which is periodic in logij^t. Since |</)(t) — (/>(u)| = 
\E{Z{t) — Z{u))\ < d-w {Z (t) , Z (u)) , (jl.27p implies that (p is continuous, and Lipschitz on 
compact intervals. □ 

Remark 1.4. As remarked above, Z(t) +log]^/^t is integer-valued. Moreover, for every 
integer k, 

P(Z(t) + logi/„ t<k)= P(Z(t) < A; - logi/„ t) = F{k - logi/„ t) 

= E{Z<k- logi/„ t) = P(Z + logi/, t < k) 
= ¥{\Z + \ogy^t\<k). 

Hence, for every t > 0, Z{t) = \Z + logij^ t \ — log^/ct t. General families of random variables 
of this type are studied in |10) . In particular, [101 Theorem 2.3] shows how (j){t) := EZ{t) in 
Theorem 11.31 can be obtained from the characteristic function of the distribution F of Z. 

Remark 1.5. The very slow convergence rate O (log~^t) in (jl.26p is because we allow 5n 
to tend to slowly. In typical applications, (5„ = for some a > 0, and then better 
convergence rates can be obtained. We have, however, not pursued this. 

Remark 1.6. Note that F and ip are influenced by the distribution of Yn for small n > 2, for 
example Is and I4; hence there is no hope for a nice explicit formula for F or depending 
only on asymptotic properties of Yn- 

Remark 1.7. The general problem of studying the number of steps until absorption at 1 
of a decreasing Markov chain on {1,2, . . . } appears in many other situations too, usually 
with quite different behaviour of Yn and X„. As examples we mention the recent papers 
studying random trees and coalescents by Drmota, Iksanov, Moehle and Roesler |3], Iksanov 
and Mohle [9j and Gnedin and Yakubovich [8]; in these papers the number killed in each 
round is much smaller than here and thus X„, is larger, of the order n or n/logn; moreover, 
after normalization Xn has a stable limit law. 
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2 Extensions 



We have assumed that we repeat the ehmination step until only one player remains. As a 
generalization we may suppose that we stop when there are at most a players left, for some 
given number a. 

Theorem 2.1. Consider the leader election algorithm described in SectionUl but stopping 
as soon as the number of remaining players is at most a, for some fixed a> 1. Suppose that 
Condition \l.l\ is satisfied. Then, the conclusions of Theorem \1.3\ hold, for some F and (p that 
depend on the threshold a. 

Proof. This generalization can be obtained from the version in Section [T] by replacing Yn by 

Suppose that Condition 11.11 holds for (Yn). It is easily seen that then Condition 11.11 holds 
for (1^) too, with the same a; for 

]E|y„ - Y^\ < aP(y„ <a) = 0{n-^~') 

and thus EY^+i — KY^ = Ey^+i — KYn + o(n^^). Consequently, Theorem 11.31 applies to {Y^), 
and the result follows. □ 

In this situation, it is also interesting to study the probability vrj(n) that the procedure 
ends with exactly i players, starting with n players; here i = 1, . . . , a and Yli=i ^i(^) = 1- 
We have a corresponding limit theorem for vrj(n). 

Theorem 2.2. Suppose that Condition holds and that a > 1 is given as in Theorem \2.1\ 
Then, 

7rj(n) = V'i(ra) + o(l), i = l,...,a, (2.1) 

for some continuous functions ipi{t) on (0, oo) which are periodic in log;^/^^, i.e. ipi{t) = 
tpi{at), and locally Lipschitz. 

Proof. A modification of the proof of Theorem 11.31 now taking Xn '■= vrj(n) and dn := 
\xn+i —Xn\ and replacing the random Xt by x\t\^ = vri([tj), yields 7rj(n + 1) — 7rj(n) = 0(l/n) 
and TTi{a~H) — TTi{a~^^~^^h) = 0{j^^^'^); hence, for any t > 0, 7ri{a~H) ijjiit) for some 
ipi{t), which easily is seen to satisfy the stated conditions. We omit the details. □ 

More generally, there is a similar result on the probability that the process passes through 
a certain state; this is interesting also for the process in Section [T] with a = 1. 

Theorem 2.3. Suppose that Condition holds and that a > 1 is given as above. Let'Ki{n), 
i >!, be the probability that, starting with n players, there exists some round with exactly i 
survivors. Then 

^i(n) = ^,(n)+o(l), i = l,2,..., (2.2) 

for some continuous functions tpiit) on (0, oo) which are periodic in logi/^t, i.e. ipi{t) = 
ipi{at), and locally Lipschitz. 



(ii) , note that Condition ll.]](iii) implies that 
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Proof. For i < a, this TTi{n) is the same as in Theorem 12.21 for each i > a, this 7rj(n) is 
the same as in Theorem 12.21 if we replace a by i. □ 

Remark 2.4. Another variation, which is natural in some problems, is to study a non- 
increasing Markov chain on {0, 1, . . . } and ask for the number of steps to reach 0; in other 
words, the time until all players are killed. In this case, we thus assume that < 1^ < n. 
This can obviously be transformed to our set-up on {1,2,...} by by increasing each integer 
by 1; in other words, we replace 1^ by 1^ := Yn-i + 1, n > 2; we can interpret this as adding 
a dummy player that never is eliminated, and continuing until only the dummy remains. 
If Condition 11.11 holds for 1^, except that 1^ = is allowed and P(yi = 0) > 0, then 
Condition 11.11 holds for too, and thus our results hold also in this case, with Xn now 
defined as the number of steps until absorption in 0. (To be precise, Xn = -^^^+1; with 
(yX'^) corresponding to {Yn)-, since we add a dummy, but there is no difference between the 
asymptotics of ^"^^ -^n-) 



3 Examples 

Example 3.1 (a toy example). For a simple example to illustrate the theorems above, let, 
for n>2,Yn= [{n + /)/2j, where I ~ Be(l/2) is or 1 with P(/ = 1) = 1/2. In other 
words, we toss a coin and let Yn be either \n/2\ or \n/2\ depending on the outcome. (If n 
is even, thus always Yn = n/2.) Note that El^ = n/2, n > 2, and that Condition 11.11 holds 
trivially, with a = 1/2. If we start with Nq = n players and m2^ < < (m -|- 1)2-', m > 1, 
then the number Nj of survivors after j rounds satisfies m < Nj < m + 1 and lEA'^- = 2^^n 
(by induction on j). Consequently, if m2^ < n < (m -|- 1)2-', 

W{Nj =m)=m + l- 2~^n, F{Nj = m + I) = 2"% - m. (3.1) 

Taking m = 1, this shows that if 2^ < n < 2^^^, then 

F{Xn = j) = = 1) = 2 - 2-%, F{Xn = i + 1) = 1 - F{Xn = j) = 2-^n - 1. (3.2) 

Hence, ([112]) holds exactly, ¥{Xn < k) = F{k - log2 n), for all A: € Z and n > 1, with 



F{x) 



0, x< -1, 

2-2-^, -l<x<0, 
[1, x>0, 



and (I1.16P holds exactly, EX„ = log2 n + </>(«), with 

(/.(2^) = 2^-W - {x-[x\) -1. 

Suppose now, as in Section [21 that we stop when there are at most a = 3 players left. 
(Similar results are easily obtained for other values of a.) If 2* < n < 2*"''^ with i > 1, we 
take j = i — 1 and note that Nj G {2, 3, 4}. If Nj = 4, the procedure ends after one further 
round with 2 players left; otherwise it ends immediately with Nj = 2 or 3 survivors. Taking 
m = 2 or m = 3 in (|3.ip . we thus find 

TT^in) = P(iVi_i = 2)+ P(iVi_i = 4) = \2^-'n - 3|, 2' < n < 2*+^ 
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Consequently, vr2(n) = ip2{n) exactly, for all n > 2, with 

^2(2^) = |2^+^-L^J -3|. 
Further, '4)-i{t) = 1 - '4)2{t) and = 0. 

Example 3.2 (a counter example). The procedure in Example 13.11 is almost deterministic. 
In contrast, the very similar but completely deterministic 1^ = [n/2j , n > 2, does not satisfy 
Condition ll.]](ii)[ In this case, Xn = Uog2 n\ and P(X„ < k) = F{k — log2 n), for all A; G Z 
and n > 1, where F{x) = l[x > 0], the distribution of Z := 0; this limit F is not continuous 
so the conclusions of Theorem 11.31 do not all hold. 



Example 3.3. A leader election algorithm studied by Prodinger [18], Fill, Mahmoud and 
Szpankowski [5] Knessl [T2], and Louchard and Prodinger [16], see also Szpankowski [ 19|, 
Section 10.5.1], is the following: Each player tosses a fair coin. If at least one player throws 
heads, then all players throwing tails are eliminated; if all players throw tails, then all survive 
until the next round. 

Except for the special rule when all throw tails, which guarantees that at least one player 
survives each round, the number Yn of survivors in a round thus has a binomial distribution 
Bi(n, ^). More precisely, if Wn is the number of heads thrown, 

Yn = Wn + nl[Wn = 0] with Wn ~ Bi(n, i). (3.3) 

Note that 

E|y„ - n/2f = E\Wn - n/2f = 0{n^), 
so (|1.4p holds for p = 6 (and, indeed, for any p > 0), and thus (|1.3p holds. Similarly, 

EYn = EWn + 2-"n = ^n + n2-". 



Thus, conditions (ii) and (iii) in Theorem 11.31 are satisfied. Also the monotonicity condition 



(i) is satisfied, because ifl<A;<n — 1 (other cases are trivial), then 



HYn+i <k)= F(l < Wn+l <k) = iP(l <Wn<k) + iP(0 <Wn<k-l) 

= <k) + \ {nWn = 0) - ¥{Wn = k)) (3.4) 

< < k). 

Hence Condition 11.11 is satisfied with a = 1/2 and Theorem 11.31 applies. 

In this case, Prodinger [18], see also Fill, Mahmoud and Szpankowski [5j, found an exact 
formula for the expectation EX„ and asymptotics of the form (jl.l6p with the explicit function 

'?k'K 

m = \ - (log 2)-^ C(l - - Xfc)e'''^'^°s^*, Xfc := r^i- (3-5) 

Fill, Mahmoud and Szpankowski [5j further found asymptotics of the distribution that can 
be written as (|1.12p with 

Fix) = ^ , (3.6) 

^ ' exp(2-^) - 1' ^ ' 
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which thus is the distribution function of Z in this case. Second and higher moments are 
considered by Louchard and Prodinger [16j. 

Prodinger [18] considered also the possibihty of stopping at a = 2 players, and showed 
(jl.l6p above in this case too, with an explicit formula for the function (of the same type 
as daSD for a = 1). 

Example 3.4. A variation of Example 13.31 studied by Janson and Szpankowski |llj . Knessl 
[12j and Louchard and Prodinger [16j is to let the coin be biased, with probability p G (0, 1) 
for heads (=survival). Then (|3.3p still holds, but with Wn ~ Bi(n,p). Conditions ll.l|nij|^iii) 
hold as above, with a = p, but arguing as in (j3.4p we see that Condition ll.]|FI)] holds if 
p > 1/2, but not for smaller p. Hence, Theorem 11.31 applies when p > 1/2. 

In fact, the results of [llj show that the conclusions ()1.12p and (I1.16P hold for all p £ (0, 1), 
for some functions F and (p explicitly given in [llj . (See also ^16j, where further higher 
moments are treated.) This suggests that Theorem 11.31 should hold more generally. Note that 
although Condition [T^i)] does not hold for p < 1/2, the difference P(y„+i < A;) - P(y„ < k) 
is at most P(l^ = 0) = (1 — p)", and thus negative or exponentially small for large n. It 
seems likely that Theorem 11.31 can be extended to such cases, by allowing a small error in 
Condition ll.]|(i)[ this then would include this leader election algorithm with a bisaed coin for 
any p E (0, 1). However, we have not pursued this. 

It was left as an open question in [llj whether for each p G (0, 1) the limit function F is 
monotone, and thus a distribution function, which means that there exists a random variable 
Z such that ()1.13p holds. By the discussion above. Theorem 11.31 shows that this holds for 
p > 1/2, but the case p < 1/2 is as far as we know still open. Cf. Remark 14.21 and Figure 
[6] below, which show that monotonicity fails in a related situation. Numerical experiments, 
based on [16, Prop. 3.1], indicate that F is monotone, at least for some choices of p < 1/2. 

A further variation of Example 13.31 is to let the probability p depend on n. The case p = 
1/n is studied by Lavault and Louchard [13j; in this case El^ is bounded and Condition ll.il 
does not hold, so Theorem 11.31 does not apply. 

Example 3.5. The special rule in Examples 13.31 and 13.41 for the exceptional case when all 
throw tails is of course necessary to prevent us from killing all players, but as we have seen, 
it complicates the analysis, especially for p < 1/2 when it destroys stochastic monotonicity 
of Yn- Note that this rule typically is invoked only towards the end of the algorithm, when 
only a few players are left. We regard the rule as an emergency exit, and it could be replaced 
by other special rules for this case. For example, an alternative would be to switch to some 
other algorithm that is fail-safe although in principle (for large n) slower; for our purpose 
this means that the present algorithm terminates, so we may describe this by letting Yn = 1 
in this case, i.e., (j3.3p is replaced by Yn = Wn + l[Wn = 0] = max(W„, 1). Note that for this 
version. Condition 11.11 holds for every p £ (0, 1), with a = p, so our theorems apply. 

An equivalent way to treat this version is to add (as in Remark 12. 4p a dummy, which 
is exempt from elimination, and to eliminate everyone else that throws a tail; we then stop 
when there are at most 2 players left (the dummy and, possibly, one real player). This is 
thus the version in Section [21 with a = 2 and y„ = 1 + Wn-i, Wm ~ Bi(m,p) (and starting 
with n + 1 players). Again, Condition 11.11 holds, with a = p, and the results in Section [2] 
apply. In particular, since invoking the special rule corresponds to eliminating everyone 
except the dummy, the probability that we have to invoke the special rule is the same as the 
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probability that the dummy version ends with only the dummy, i.e. 7ri(n+ 1) in the notation 
of Theorem \2.2\ the asymptotics of this probability is thus given by (12. ip . 



Example 3.6. Prodinger [T7] (for p = 1/2) and Louchard and Prodinger [15] studied a 
version of Examples 13.31 and 13.41 where, as in Remark 12.41 we allow all players to be killed 
and let X„ be the time until that happens. Thus, in each round, each player tosses a coin 
and is killed with probability 1—p (and we do not have any special rule). Additionally, there 
is a demon, who in each round kills one of the survivors (if any) with probability v £ [0, 1]. 
We thus have the modification discussed in Remark 12^ with 



Yn = max{Wn -1,0), Wn ^ Bi(n,p), ~ 

where Wn and ly are independent; thus also = Yn-i + 1 = max(W„_i + 1 — 1^,1). 
Condition 11.11 holds, in the modification for absorption in 0, and thus our results apply. In 
fact, Louchard and Prodinger [T5| show (I1.12P and (I1.16P for this problem with explicitly 
given F and (p; they further give an extension of (jl.l6p to higher moments. 

As remarked in [iTl US], the special case v = 1 \s equivalent to approximate counting 
and V = \s equivalent to the cost of an unsuccessful search in a trie; in the latter case, 
Xn is simply the maximum of n i.i.d geometric random variables which can be treated by 
elementary methods. 

Example 3.7. W.R. Franklin [7] proposed a leader election algorithm where the n players 
are arranged in a ring. Each player gets a random number; these are i.i.d. and, say, uniform 
on [0, 1]. (Since only the order of these numbers will matter, any continuous distribution will 
do; moreover, it is equivalent to let . . . , ^„ be a random permutation of 1, ... , n.) A player 
survives the first round if her random number is a peak; in other words, if ^i, . . . , are i.i.d. 
random numbers, then player i survives if > and ^ > ^j+i (with indices taken modulo 
n). We may ignore the possibility that two numbers and are equal; hence we may as 
well require > ^i-i and > ^j+i- 

In Franklin's algorithm, the survivors continue by comparing their original numbers in the 
same way with the nearest surviving players; this is repeated until a single winner remains. 
We have so far not been able to analyse this algorithm. It is easy to verify that even if 
we condition on the number m of survivors after the first round, the ml possible different 
orderings of the survivors do not appear with equal probabilities, which means that the 
algorithm is not of the recursive type studied in this paper. (For example, starting with 
a ring of 8 players and conditioning on having 4 survivors (peaks) in the first round; the 
probability of getting 2 survivors in the second round is 10/34, and not 1/3 as in the uniform 
case.) 

However, we can study a variation of Franklin's algorithm, where the survivors draw 
new random numbers in each round. This is an algorithm of the type studied in Section [H 
with Yn given by the number of peaks in a random permutation, regarded as a circular 
list. Note that there is always at least one peak (the maximum will always do), so y„ > 1 
as required. It is easily seen that inserting a new player will never decrease the number 
of peaks; hence Yn is stochastically increasing in n. Further, we have Yn = ^^h, where 
li := > max(^j_i, ^j+i)] is the indicator that player i survives (again, indices are taken 



modulo n). If n > 3, then E/j = 1/3 by symmetry, and thus El^ = n/3. In particular, (ii) 



holds with a = 1/3. Furthermore, Ii and Ij are independent unless |i — j| < 2 (mod n), and 
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similarly for sets of the indicators /j, and it follows easily that — E^)^ = 0{n^), and 
thus (jl.4[) holds for p = 6. (Indeed, (11. 4p holds for all even p by this argument.) Consequently, 
Condition 11.11 holds and Theorem 11.31 applies, with 1/a = 3. 

Remark 3.8. It might be shown that for the true Franklin algorithm, the expected number 
of survivors after 2 rounds is C2n + o(n), where 

C2 = ~ ^f^f + « 0.1096868681. (3.7) 
384 

(Details might appear elsewhere.) In comparison, for the variation with new random numbers 
each round, it is easily seen that the expected number after k rounds is (l/3)'^'n + o(l), for 
any fixed k; in particular, after two rounds it is n/9 + o{n). Note that C2 in (|3.7p is slightly 
smaller than 1/9, and thus better in terms of performance. It might have been hoped that 
the true Franklin algorithm is asymptotically equivalent to the variation studied here, but 
the fact that C2 ^ 1/9 suggests that this is not the case. Nevertheless, we conjecture that 
Theorem 11.31 remains true for the Franklin algorithm, for some unknown a < 1/3. 

Example 3.9. In both versions in Example 13.71 the players are arranged in a circle. Al- 
ternatively, the players may be arranged in a line. We use the same rules as above, but we 
have to specify when a player at the end (with only one neighbour) is a peak. There are two 
obvious possibilities: 

(i) Never regard the first and last players as peaks. (Define = Cn+i = +oo.) 

(ii) Regard them as peaks if .^i > ^2 and > S,n-i, respectively. (Define .^0 = S,n+i = —00.) 

In the first case, it is possible that there are no peaks, and thus we have to add an emergency 
exit as in Example 13.51 

As in Example 13.71 there are two versions (for each of (i) and (ii)): we may use the same 
random numbers in all rounds, or we may draw new ones each round. 

In the latter case, we are again in the situation of Section [TJ In both cases (i) and (ii), the 
distribution of Yn is related to the distribution in the circular case in Example 13.71 Indeed, 
if we start with a circular list of n + 1 numbers and eliminate the player with the largest 
number, then the remaining n numbers form a linear list, and the peaks in this list using 
version (i) equal the peaks except the maximum one in the original circular list. Similarly, if 
we instead eliminate the player with the smallest number, then the peaks in the remaining 
list using version (ii) equal the peaks in the original list. Hence, if Zn is the number of peaks 

in a random circular list of length n, then (i) yields Yn = Zn+i — 1 and (ii) yields Yn = Zn+i- 
In both cases, this implies that Condition 11.11 holds (provided we add a suitable emergency 
exit in case (i)) because it holds for Z„. Consequently, Theorem 11.31 applies to both these 
linear versions of (the variation of) Franklin's algorithm, again with 1/a = 3. 

4 First variation of the Franklin leader election algorithm. 
The linear case 

We assume that the survivors draw new random numbers in each round and that they are 
arranged in line. We use possibility (i) of Example 13.91 We start with a set of n players. We 
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assign a classical permutation of {1, . . . ,n} to the set, all players corresponding to a peak 
stay alive, the other ones are killed. If there are no peaks, we choose the following emergency 
exit: a player is chosen at random (this is assumed to have cost), indeed in the original 
game, one deals with circular permutations, so there always exists at least one peak, here we 
approach the problem with a classical inline permutation. 

What is the distribution of the number Xn of phases (or rounds) before getting only one 
player? 



4.1 The analysis 

Let 



Yn 

P{n,k) 

n(n,j) 



number of peaks, starting with n players, 
¥[Yn = k] = ¥[k peaks, starting with n players], 

P(X„ = j) = P[j phases are necessary to end the game, starting with n players], 

j 

A(n, j) := n(n, k) = P[at most j phases are necessary, starting with n players]. 

A:=0 



We will sometimes use the subscript 1 to distinguish these from the circular case discussed in 
Section [5l 

First of all, we know (Carlitz [2]; see also (BJ Chapter 3]), that the pentavariate generating 
function (GF) of valleys (uq), double rise (wi), double fall (u'l) and peaks (^2) is given by 

^, . 6 vi + 5tan(z5) vi 

I{z,u) = —- . , 

U2 — vi tan[zdj U2 

with 



vi = (ni + u[)/2, 6 = ^JuQU2 - vf. 
This gives the GF of the number of peaks: 

tan[z(M - 1)1/2] 



(4.1) 



(u- 1)1/2 _tan[z(n- 1)1/2]' 

hence the mean M and variance V of the number of peaks, for n > 2 and n > 4, respectively: 

M(n) = (n - 2) /3, V(n) = 2(n + 1) /45. (4.2) 

This GF is also given in Carlitz [2]. Moreover, from [6l Chapter 9], we know that the 
distribution P is asymptotically Gaussian. This is also proved in Esseen [3] by probabilistic 
methods. 

Let x{n) be the mean number of phases, E(X„), starting with n players. As we shall see, 
the initial values are 

j;(0) = x(l) = 0, x(2) = j;(3) = x(4) = 1. 

Since (14. 2p yields M{n) + 1 = (n + l)/3 for n > 2, we have (approximating by using this for 
n < 1 too) that the mean number of players c(j) still alive after j phases is 

c(j)~3-^(n + l)-l. 
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Table 1: Pi(n, k) 

(An induction easily yields the exact formula |c(j) — 3~-^n| < 1 for all j.) If we want c(j) = 1, 
this leads to the approximation 

x{n) K, j K, loggn - log3 2. 

We see from Theorem 11.31 (which applies by Example I3.9p that this is roughly correct, but 
the constant — log3 2 has to be replaced by a periodic function (/>(n). 
Let us now construct 11. We have, for n > 2 and j > 1, 

L(n-1)/2J 

n(n,j)= P{n,m{k,j-l). (4.3) 

fc=0 

We have the initial values 

n(o,o) = i, n(o,j) = o,j >o, n(i,o) = i, n(i,j) = o,i > o, n(n,o) = o,n > 2. 

Also 

n(2,i)=n(3,i) = n(4,i) = i. 

Some values of P and 11 is given in Tables [1] and [2j (We use in the tables and figures the 
subscript 1 to emphasise that we deal with the linear case.) 
Denoting the jih. column of 11 by ■k'^^\ we have 

For j > 2, it suffices to consider /c > 2 in (j4.3p . so we need only the matrix {P{n,k))n^k>2- 
Since P{n, k) = ii k > {n — l)/2, this matrix is triangular, and so is P^~^ . 

But 7r(^)(n) < 10^'^, n > 20, so numerically, the significant columns of P^^^ are the first 
20 columns. Also, we see the importance of the initial first column of 11. Moreover, for 
n > 75, P{n, k) is indistinguishable from the Gaussian limit. So we have used the expansion 
of the GF (14. ip for n < 75 and the Gaussian limit afterwards in our numerical calculations. 
Of course we have 

oo 

x{n) = ^n(n, j)j, 

j=0 

and 

oo 

x{n) = l + YP{n,k)x{k), n > 2. 

fc=2 
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Table 2: Il\{n,j) 

Remark 4.1. Another approach could be the following: Let 

I{k) := \ one of the phases has k players |, 
/(j, k) := I phase j has k players |, 

Qik) ■.= F[I{k) = l], 
R{j,k) :=F[I{j,k) = l], 

n—k 

Q{k)=¥[\/I{j,k) = l\ =Y,Rij,k), s.s¥[I{j,k)M{i,k)]=0, if i^j, 
j i=i 

RU, k) = Y, RU - 1, l)Pil. k), j > 1, and i?(0, k) = dkn 
I 

n n 

x{n) =W.^I{k)\ = Y,Q{k). 

1 1 

A plot of x{n) — log3 n versus log3 n is given in Figure [5] for n = 50, . . . , 500. The 
oscillations expected from (jl.l6p are clear. 

Recall that according to Theorem II. 3 [ there exists a limiting distribution function F[x) = 
F\{x) (in a certain sense) for X^- In Figure[3l we approximate this distribution function F{x) 
by plotting A(n, j) = P(X„ < j) against j — logs n for n = 20, . . . , 500, cf. (I1.12p . We have 
also plotted a scaled Gumbel distribution; the fit is bad. 

Similarly, in Figure [4] we show the probability n(n, j), n = 150, . . . , 500, plotted against 
j — log3 n. The fit with a Gaussian distribution is equally bad. 

The few scattered points of both figures are actually due to small n and the propagation 
of the more erratic behaviour for n = 1, . . . , 40 shown in Figure O 

So we observe the following facts: 

(i) A first regime (n = 1, . . . ,40) create some scattered points which almost look like two 
distributions. 

(ii) Between n = 40 and n ~ 75 , a limiting distribution is attained; n(n, j) is concentrated 
on j = log3n + 0(1), 
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Figure 2: x\{n) — log^ n versus log3 n, n = 50, . . . , 500 



(iii) For n > 75, the limiting Gaussian for P with its narrow dispersion, intuitively 
induces, with (j4.3p . a propagation, with some smoothing, of the previous distribution. 
We attain the limiting distribution F{x) given by Theorem 11.31 

(iv) At most two values carry the main part of the probability mass Il{n,j). This is clear 
from the observed range of F{x) in Figures [3] and HI 

(v) P is triangular, and so is P^ . Also, P^{n,k) = 0(1), with k = 0(1), only if j = 
log3n + 0(l). 

(vi) As F{x) is absolutely continuous, we can derive, as in [T3], modulo some uniform 
integrability conditions, all (periodic) moments of in particular x{n). 

(vii) The effect of initial values is now clear. To illustrate this we have changed to 11(0, 1) = 
n(l,0) = 1, which means that we add a cost 1 for the extra selection required when 
the algorithm terminates with no element left. This leads to Table [3l The equivalent of 
Figures El [Hand [5] is given in Figures [HIT] and El Note that Xn no longer is stochastically 
monotone in n; we have by definition Xq = 1 > = Xi, and Table El shows other 
examples of non-monotonicity for small n. Moreover, Figure El shows that the non- 
monotonicity persists for large n; we clearly have convergence to a limit function, G{x) 
say, but the limit is not monotone and thus not a distribution function as in Figure HI 
and, more generally, in Theorem 11.31 



Remark 4.2. Note that the example in (vii) and Figure El does not contradict Theorem 11.31 



because Xn now is defined with other initial values than in Theorem 11.31 Nevertheless, it is 
a warning that monotonicity of the limit should not be taken for granted in cases such as 
Example 13 .41 with p < 1/2 where the monotonicity assumption of Theorem ll.3l is not satisfied. 
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o : observed 

— : Gumbel distribution 
Figure 3: Ai(n, j) = P(X„ < j) versus j — log^, n, approximating Fi{x), n = 20, . . . , 500 
4.2 Periodicities 

Let ip{a) := / e'^^ f{x) dx be the Laplace transform of the limiting distribution F. (We do not 
know whether ■(/'(o^) exists in general, although we conjecture so, but we really only need it for 
imaginary a, i.e., the characteristic function of F.) Similarly, let 1^(0) := J e°^^ AF{x) dx be 
the Laplace transform of AF. Since AF(x) := F{x)—F{x—1) = f{x—t) dt, AF = f *X[o^i], 
and thus, since X[o,i] Laplace transform (e" — l)/a, 

- 1 

a 

With the usual machinery (see [13] or [10\ Theorem 2.3]), we obtain from Theorem 11.31 and 
Remark [131 assuming some technical conditions that are very likely to hold but not rigorously 
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Table 3: n(n,j), other initialization 
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Figure 4: ni(n, j) = P(X„ = j) versus j — log3 n, approximating AFi(x), n = 150, . . . , 500 
verified, 

x{n) := KXn = log^ n + rhi + t(;i(log3 n) + o(l), (4-4) 

where 

rhi := ^'(0) = V'(0) + i, (4.5) 
u;i(x) := y V;'(27r/i)e2-'- = V ^l^g^-'-. (4.6) 

With the observed values of n(n, j) (see Figure we have computed the Laplace trans- 
form numerically (with a variable step Euler-MacLaurin) as follows. Assume that we have N 
computed values {Il{ni,ji) : i = 1, . . . ,N}. Setting x = j — logn, this gives {n(nj,Xj) : i = 
1, . . . , N}. Sorting wrt Xj, we write this as {n(nfc, y^) : Uk < Vk+i, k = 1, . . . , N}. Construct 
a numerical Laplace transform 

i^ia) = ^e"2"=n(nfc,yfc)(yfc+i -yfe_i)/2. 

Using this numerically computed ip in ()4.4p - (j4.6p . we compute rhi + 'Wi{x), which fits quite 
well with the observed periodicities of x{n) — log3n in Figure [21 the comparison is given in 
Figure [9l 



20 





- 


-1 

0.3 








0.6 








-0.4 








0.2 




-1.5 


1 -0.5 





0.5 



Figure 5: Ili{n,j) versus j — log^ n, n = 1, . . . , 40 

5 Second variation of the Franklin leader election algorithm. 
The circular case 

If we denote by Pc{n,k) the distribution of the number of peaks in the circular case and 
by P\{n,k) the distribution in the hnear case, we know, by Example 13.91 that Pc{n.,k) = 
Pi{n — l,k — 1). It is easy to check that this leads to, for n > 3 and n > 5, respectively, 

M(n) = n/3, V(n) = 2n/45, 

which also is easy to see probabilistically, by writing Yn as the sum of the n indicators 
[player z is a peak|, and noting that indicators with distance at least 3 are independent.) 
The initial values are now given by 

1) = 1, p(2, 1) = 1, p(3, 1) = 1, n(o, 0) = 1, n(i, 0) = 1, n(2, i) = i, n(3, i) = i. 

The corresponding pictures are given in Tables [H and \5\ (We use a subscript c for the circular 
case.) 

A plot of Ac(n,j) versus j — log3n for n = 20, . . . ,500 is given in Figure [lOl We see 
that there are fewer scattered points than in the linear case. Let us mention that the fits 
with Gumbel or Gaussian are equally bad. A comparison of Ac{n,j) with Ai(n,j) is given in 
Figure [m No numerical relation exists between the two distributions. 

nc(ra, j) versus j — log3n, n = 1, . . . ,40, is plotted in Figure [T2l These initial points are 
now scattered less than in the linear case. 

The observed versus computed periodicities are given in Figure [T3l 

In conclusion, apart from numerical differences, the behaviour of our two variations are 
quite similar. 

Note that the mean number of needed messages is asymptotically 2nlog3(n), as we use 
2n messages per round. Franklin gives an upper bound 2nlog2(?T.). 
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Table 5: Uc{n,j) 
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Figure 6: A(n, j) versus j — log^ n, other initialization, n = 5, . . . , 500 
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Figure 8: Il(n,j) versus j — log^ n, n = 1, . . . , 100, other initialization 
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Figure 9: observed x{n) — log3 n (o) and computed with ()4.4I) (line) periodicities versus log3 n 
(linear case), n = 50, . . . , 500 
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Figure 10: Ac(n, j) versus j — log^ n, approximating Fc{x), n = 5, . . . , 500 
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Figure 11: A comparison of Ac(n, j) (right curve) with Ai(n, j) (left curve) 
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Figure 12: Ilc{n,j) versus j 
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